High real-space resolution measurement of the local structure of Gai_ 

x-ray diffraction 



c In x As using 



On 

In 



o 

CZ3 



Ctf 



c 

O 

o 



> 

On 
On 
O 
^D 
O 
On 
On 

-i— > 
03 



i 

-a 

c 
o 
o 



s 



V. Petkov 1 , I-K. Jeong 1 , J. S. Chung 1 , M. F. Thorpe 1 , S. Kycia 2 and S. J. L. Billinge 1 

1 Department of Physics and Astronomy and Center for Fundamental Materials Research, Michigan State University, East 

Lansing, MI 48824-1116. 2 Cornell High Energy Synchrotron Source, Cornell University, Ithaca, NY 14853 

(February 7, 2008) 

High real-space resolution atomic pair distribution functions (PDF)s from the alloy series 
Gai-aJnzAs have been obtained using high-energy x-ray diffraction. The first peak in the PDF 
is resolved as a doublet due to the presence of two nearest neighbor bond lengths, Ga-As and In- As, 
as previously observed using XAFS. The widths of nearest, and higher, neighbor pairs are analyzed 
by separating the strain broadening from the thermal motion. The strain broadening is five times 
larger for distant atomic neighbors as compared to nearest neighbors. The results are in agreement 
with model calculations. 



The average atomic arrangement of crystalline semi- 
conductor alloys is usually obtained from the position 
and intensities of the Bragg peaks in a diffraction ex- 
periment pi, and the actual nearest neighbor and some- 
times next nearest neighbor distances for various pairs of 
atoms by XAFS measurements M . In this Letter we show 
how high energy x-ray diffraction and the resulting high- 
resolution atomic pair distribution functions (PDF)s can 
be used for studying the internal strain in Gai-^In^As 
alloys. We show that the first peak in the PDFs can 
be resolved as a doublet and, hence, the mean position 
and also the widths of the Ga-As and In- As bond length 
distributions determined. The detailed structure in the 
PDF can be followed out to very large distances and the 
widths of the various peaks obtained. We use the con- 
centration dependence of the peak widths to separate the 
strain broadening from the thermal broadening. At large 
distances the strain broadening is shown to be about five 
times larger than for nearest neighbor pairs. Using a sim- 
ple valence force field model, we get good agreement with 
the experimental results. 

Ternary semiconductor alloys, in particu- 
lar Gai-^In^As, have technological significance because 
they allow important properties, such as band-gaps, to be 
tuned continuously between the two end-points by vary- 
ing the composition x. Surprisingly, there is no com- 
plete experimental determination of the microscopically 
strained structure of these alloys. On average, both GaAs 
and InAs form in the zinc-blende structure where Ga 
or In and As atoms occupy two inter-penetrating face- 
centered-cubic lattices and are tetrahedrally coordinated 
to each other M. However, both extended x-ray absorp- 
tion fine structure (XAFS) experiments g and theory ||] 
have shown that Ga-As and In- As bonds do not take some 
average value but remain close to their natural lengths 
of L° Gsl _ As = 2.437 A and Lf n _ M = 2.610 A in the al- 
loy. Due to the two considerably different bond lengths 
present, the zinc-blende structure of Gai_ x In x As alloys 
becomes locally distorted. A number of authors (^-||] 
have proposed distorted local structures but there has 
been limited experimental data available to date. The 
fully distorted structure is a prerequisite as an input for 



accurate band structure and phonon dispersion calcula- 
tions ||. 

The technique of choice for studying the local structure 
of semiconductor alloys has been XAFS [|]||. However, 
XAFS provides information only about the immediate 
atomic ordering (first and sometimes second coordina- 
tion shells) and all longer-ranged structural features re- 
main hidden. To remedy this shortcoming we have taken 
the alternative experimental approach of obtaining high- 
resolution PDFs of these alloys from high energy x-ray 
diffraction data. 

The PDF is the instantaneous atomic density-density 
correlation function which describes the local arrange- 
ment of atoms in a material H. It is the sine Fourier 
transform of the experimentally observable total struc- 
ture function obtained from powder diffraction measure- 
ments. PDF analysis yields the real local structure 
whereas an analysis of the Bragg scattering alone yields 
the average crystal structure. Determining the PDF has 
been the approach of choice for characterizing glasses, liq- 
uids and amorphous materials for a long time S. How- 
ever, its widespread application to study crystalline ma- 
terials has been relatively recent M . Very high real-space 
resolution is required to differentiate the distinct Ga-As 
and In-As bond lengths present in Gai-^In^As. High 
real-space resolution is obtained by measuring the struc- 
ture function, S(Q) (Q is the amplitude of the wave vec- 
tor), to a very high value of Q (Q max > 40 A -1 ). An 
indium neutron absorption resonance rules out neutron 
measurements in the Gai_ x In x As system. We therefore 
carried out x-ray powder diffraction measurements. To 
access Q values in the vicinity of 40-50 A -1 it is necessary 
to use x-rays with energies > 50 keV. The experiments 
were carried out at the A2 56 pole wiggler beamline at 
Cornell High Energy Synchrotron Source (CHESS) which 
is capable of delivering intense x-rays of energy 60 keV. 
Six powder samples of Gai-xIn^As, with x = 0.0, 0.17, 
0.5, 0.67, 0.83 and 1.0, were measured. The samples were 
made by standard methods and the details of the sample 
preparation will be reported elsewhere J10[ . All measure- 
ments were done in symmetric transmission geometry at 
10K. Low temperature was used to minimize thermal vi- 
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FIG. 1. The reduced structure factor F(Q) = Q[S(Q) - 1] 
for Gai-^In^As measured at 10K for various concentrations. 
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FIG. 2. The reduced PDF, G(r), for Gai-^In^As measured 
at 10 K for various concentrations. 



bration in the samples, and hence to increase the sensi- 
tivity to atomic static displacement amplitudes. A dou- 
ble crystal Si(lll) monochromator was used. Scattered 
radiation was collected with an intrinsic germanium de- 
tector connected to a multi-channel analyzer. The elastic 
component was separated from the incoherent Compton 
scattering before data analysis [fL0|. Several diffraction 
runs were conducted with each sample and the result- 
ing spectra averaged to improve the statistical accuracy. 
The data were normalized for flux, corrected for back- 
ground scattering, detector deadtime and absorption and 
divided by the average form factor to obtain the total 
structure factor, S(Q) |7]||[flJ. Details of the data pro- 
cessing are described elsewhere Effl. Correction proce- 
dures were done using the program RAD [ |12| . Experi- 
mental reduced structure factors, F(Q) — Q[S(Q) — 1], 
are shown in Fig. [I]. The corresponding reduced PDFs, 
G(r), obtained through a Fourier transform 
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are shown in Fig. 0. The data for the Fourier transform 
were terminated at Q m ax = 45 A -1 beyond which the 
signal to noise ratio became unfavorable. This is a very 
high momentum transfer for x-ray diffraction measure- 
ments; for comparison, Q ma x from a Cu K Q x-ray tube 
which is less than 8 A -1 . 

Significant Bragg scattering (well-defined peaks) are 
immediately evident in Fig. 1 up to Q ~ 40 A -1 in the 
end-members, GaAs and InAs. This implies that the 
samples have long range order and that there is little po- 
sitional disorder (dynamic or static) on the atomic scale. 
The Bragg-peaks disappear at much lower Q- values in the 
alloy data: the samples are still long-range ordered but 



they have significant local positional disorder. At high-Q 
values, oscillating diffuse scattering is evident. This has 
a period of 27r/2.5 A -1 and contains information about 
the shortest atomic distances in Gai-^In^As alloys seen 
as a sharp first peak in G{r) at 2.5 A (see Fig. g). In 
the alloys, this peak is split into a doublet as is clearly 
evident in Fig. [2| with a shorter Ga-As bond and a longer 
In- As bond. This peak is shown on an expanded scale in 
the inset to Fig. |] for all the compositions measured. We 
determined the positions of the two subcomponents of 
the first PDF peak, i.e. the mean Ga-As and In- As bond 
lengths, and the results are shown in Fig. ||. Also shown 
is the room temperature result previously obtained in 
the XAFS study of Mikkelson and Boyce f|. There is 
clearly good agreement. The PDF-based bond lengths 
are shifted to smaller lengths by about 0.012 A since our 
data were measured at 10K, whereas the XAFS experi- 
ments were at room temperature. The nearest neighbor 
peak is the only peak which is sharp in the experimental 
PDFs as can be seen in Fig. 2. From the second-neighbor 
onwards the significant strain in the alloy samples results 
in broad atom-pair distributions without any resolvable 
splitting. Model calculations show that this broadening 
is intrinsic and not due to any experimental limitations. 
The strain in Gai-^In^As was quantified by fitting the 
individual peaks in experimental PDFs. We used Gaus- 
sians convoluted with Sine functions which account for 
the experimental resolution coming from the finite Q ma x- 
The FWHM of the resolution function is 0.086 A. This is 
significant for the near-neighbor peaks as shown in Fig. 3, 
but is much smaller than the width of the high-r peaks. 
The high-r peaks are fit using the PDFFIT modeling pro- 
gram fl3fl assuming the virtual crystal zinc-blende struc- 
ture and refining displacement parameters. The resulting 
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FIG. 3. Solid symbols: Ga-As and In-As bond-lengths vs. 
composition as extracted from the present PDFs. Open sym- 
bols: room-temperature XAFS results from Ref. H. Inset: 
Split nearest neighbor PDF peak from the data (symbols) 
and the model (solid lines). 



mean-square Gaussian standard deviations are shown in 
Fig. 4. One can see Ga-As and In-As bond lengths are 
sharply peaked about the mean value whereas the static 
strain on more distant neighbors is five times larger than 
on the near-neighbors. The strain peaks at a composi- 
tion cc=0.5 and affects the common (As) more than the 
mixed (metal) sublattice. 

In order to better understand these results, we have 
modeled to the static and thermal disorder in the alloy 
by using a Kirkwood potential |L4j]. The key element 
in this potential is the central force term that connects 
nearest neighbor atoms and would like to keep each bond 
at its natural (unstrained) length. The potential contains 
nearest neighbor bond stretching force constants a and 
force constants (3 that couple to the change in the an- 
gle between adjacent nearest neighbor bonds. We choose 
these parameters to fit the end members |3| with aca-As 

= 96N/m, ai n _As = 97N/m, /3 G a-As-Ga = /^As-Ga-As 

= lON/m and An-As-in = /^As-in-As = 6N/m. The ad- 
ditional angular force constant required in the alloy are 
taken to be the geometrical mean, so that /3ca-As-in = 
\/(/3Ga-As-Ga-/?in-As-in)- We have constructed a series 
of cubic 512 atom periodic supercells in which the Ga 
and In atoms are distributed randomly according to the 
composition x. The system is relaxed using the Kirk- 
wood potential to find the displacements from the vir- 
tual crystal positions. The volume of the supercell is 
also adjusted to find the minimum energy. Using this 
strained static structure, a dynamical matrix has been 
constructed and the eigenvalues and eigenvectors found 
numerically. From this the Debye- Waller factors for all 



the individual atoms in the supercell can be found and 
hence the PDF of the model by including the Gaus- 
sian broadening of all the subpeaks. We have shown 
previously |15| that this is the correct procedure within 
the harmonic approximation. The model-PDF is plot- 
ted with the data in the inset to Fig. 2 and in Fig. 5. 
The agreement at higher-r is comparable to that in the 
r-range shown. All the peaks shown in the Figures con- 
sist of many Gaussian subpeaks. The overall fit to the 
experimental G(r) is excellent and the small discrepan- 
cies in Fig. 5 between theory and experiment are proba- 
bly due to small residual experimental errors. Note that 
in comparing with experiment, the theoretical PDF has 
been convoluted with a sine function to incorporate the 
truncation of the experimental data at Q m ax — 45 A. 
The technique discussed above could be extended using 
a better force constant model with more parameters, but 
does not seem necessary at this time. 

The thermal and strain contributions to the widths of 
the individual peaks in the reduced PDF act indepen- 
dently, as expected and as confirmed by our supercell 
calculations described in the previous paragraph. We 
therefore expect the squared width A to be a sum of the 
two parts. The thermal part a is almost independent of 
the concentration and we fit a 2 by a linear function of 
the composition x between the two end points in Fig. 4. 
To better understand the strain model it is convenient 
to assume that all the force constants are the same and 
independent of chemical species. Then it can be shown 
|16| for any such model that 
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where the subscripts ij refer to the two atoms that lead 
to a given peak in the reduced PDF. For the Kirkwood 
model the A^ are functions of the ratio of force constants 
/3/q only. It further turns out that the Aij are indepen- 
dent of whether a site in one sublattice is Ga or In, so 
we will just refer to that as the metal site. Taking mean 
values from the force constants used in the simulation we 
find that j3/a = 0.83, and that for nearest neighbor pairs 
A^ = 0.0712. For more distant pairs the motion of the 
two atoms becomes incoherent so that Aij = Ai + Aj 
and we find that for the metal site Ai = 0.375 and for 
the As site Aj = 1.134. The validity of the approxi- 
mation of using mean values for the force constants was 
shown to be accurate by calculating the model-PDF for 
all compositions as described above and comparing to 
the prediction of Eq. (2) ]lq] . Equation 2 shows good 
agreement with the data for near and far neighbor PDF 
peaks, and for the different sublattices, over the whole 
alloy range, as shown in Fig. 4, using only parameters 
taken from fits to the end-members. There is a consid- 
erably larger width associated with the As-As peak in 
Fig. 4 when compared to the Me-Me peak, because the 
As atom is surrounded by four metal cations, providing 
five distinct first-neighbor environments pip). The theo- 
retical curve in the lower panel of Fig. 4 is predicted to 
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FIG. 4. square of the PDF peak widths for far neighbors 
(top panel) and nearest neighbors (lower panel) separated by 
sub-lattice type. Symbols: values from the data. In the lower 
panel the open symbols are for the Ga-As and the closed 
symbols the In-As bond. See text for details. Solid lines: 
theory. The mean-square static and thermal distortions are 
added. Here Me represents both the metals Ga and In, which 
behave in the same way. Note the scale in the lower panel is 
expanded by a factor 10 compared to the upper panel. 
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FIG. 5. Experimental (open circles) and theoretical (solid 
line) PDFs for Gai-^In^As for concentrations x = 0.5 and 
x = 033. 



be the same for the Ga-As and In-As bond length dis- 
tribution, using the simplified approach. The Kirkwood 
model seems adequate to describe the experimental data 
at this time, although further refinement of the error bars 
may require the use of a better potential containing more 
parameters. 

In summary, we report for the first time a high-real- 



space-resolution measurement of the PDF of Gai_;rIn x As 
(0 < x < 1) alloys. The PDF allows the local distortions 
away from the average structure over a wide range of r to 
be determined in disordered crystals such as these. The 
nearest-neighbor Ga-As and In-As bond lengths in the 
alloys are clearly resolved. Significantly greater disorder 
exists in the more distant neighbor length distributions 
in the alloys. The experimental results are well fit over 
a wide range or r using a Kirkwood model. Because 
the agreement between theory and experiment is good at 
both short and large distances, the Kirkwood model can 
be used with some confidence to generate strained alloy 
structures for use in the calculation of electronic band 
structure and phonon dispersion curves. 
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